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摘要 :研究 生态 系统 服务 之 间 以 及 与 土地 利用 类 型 之 间 的 复杂 关系 ,进行 生态 系统 服务 导向 下 的 主 地 利用 优化 /对 实现 社会 经 
济 与 生态 保护 协调 可 持续 发 展 具有 重要 的 意义 。 以 渭河 流域 关中 -天 水 经 济 区 段 ( 简称 关 天 段 ) 为 研究 对 象 ;在 子 流域 尺度 估算 
了 2000—2013 年 的 生物 多 样 性 、 固 矶 和 产 水 3 种 生态 系统 服务 量 , 定 量化 研究 彼此 之 间 的 权衡 和 协同 关系 ;并 利用 CA-Markov 
模型 预测 2050 年 土地 利用 变化 情景 ,以 此 测算 未 来 情景 下 的 生态 系统 服务 并 引入 生产 可 能 性 边界 (PPF ) 方 法 ,进行 权衡 与 协 
同 分析 ,探究 2050 年 最 优生 态 系统 服务 下 的 土地 利用 格局 。 研 究 得 出 :2000 一 2013 年 研究 区 的 总 产 水 量 呈 增加 趋势 。 秦 岭 
与 六 盘山 地 区 总 固 碳 量 较 大 ,生物 多 样 性 质量 较 高 ;而 关中 平原 地 区 人 口 和 城市 分 布 较 多 ; 固 碳 能 力 较 差 , 生 境 退 化 度 普遍 较 
高 。@) 固 碳 与 生物 多 样 性 之 间 呈 现 协同 关系 ,而 产 水 与 生物 多 样 性 和 产 水 与 固 碳 之 间 鬼 表现 为 权衡 关系 。@) 绘 制 了 2050 年 的 
生态 系统 服务 下 的 最 优 帕 累 托 效率 曲线 ,并 进一步 进行 权衡 协同 分 析 , 得 到 相对 应 的 主 地 利用 优化 图 。 将 生态 系统 服务 及 权衡 
需求 纳入 到 土地 决策 和 优化 管理 过 程 中 ,将 有 助 于 研究 区 的 生态 ,社会 和 经 济 达到 共 赢 状态 。 
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Trade-offs and optimization among ecosystem services in the Weihe River basin 
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Abstract: With the rapid increase in human population and the excessive utilization of natural resources, humans frequently 
enhance the production of some services at the expense of others. Therefore, it is important to accurately understand the 
trade-offs and synergies among ecosystem services and explore land use optimization in order to improve ecosystem services 
in the Weihe River watershed. In this study, we first simulated and mapped three key ecosystem services at the sub-basin 
level from 2000 to 2013, including water yield, carbon sequestration, and biological diversity. Then, we quantitatively 
analyzed the dynamic changes in trade-offs and synergies between ecosystem services by using R language programs and 
production possibility frontiers ( PPF). Finally, the CA-Markov model was used to predict land use change scenarios for 
2050. We estimated and optimized the qualities of various ecosystem services in 2050 and the optimal pattern of land use 
was generated based on the Optimal Pareto Frontier between ecosystem services. The results indicate that; () the water yield 
increased in the Weihe River watershed from 2000 to 2013. The total quantity of carbon sequestered in the Qinling 
Mountains and Liupan Mountains was relatively large and the biodiversity quality was relatively high. The Guanzhong Plain 
region had a large population and urban distribution, the carbon sequestration capacity was poor, and the habitat 


degradation degree was generally higher. @There were tradeoffs between water yield and carbon sequestration, as well as 
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water yield and biological diversity. There was synergy between carbon sequestration and biological diversity. @) The optimal 
pattern of land use was generated based on the Optimal Pareto Frontier between ecosystem services in 2050. Our research 
suggested that integrating ecosystem services and trade-off analyses into land system architecture and optimal management is 


of great value to the sustainable development of the environment and economy. 
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生态 系统 为 人 类 提供 赖 以 生存 的 目 然 环境 条 件 及 效用 ,人 类 直接 或 间接 地 从 中 获取 所 有 惠 益 就 是 生态 系 
SURE" 。 生 态 系统 服务 种 类 多 种 多 样 , 不 仅 可 以 为 人 类 的 生存 直接 提供 粮食 、 氧 气 `, 水 等 ,并 具有 净化 污 
Qe 涵养 水 源 .保护 生物 多 样 性 水土 保 持 等 功能 ”。 生 态 系统 服务 之 间 关 系 复杂 ,相互 影响 ;具体 表现 为 此 
消 披 长 的 权衡 关系 和 相互 增益 的 协同 关系 。 正 确 理解 生态 系统 服务 间 的 权衡 与 协同 关系 对 提高 人 类 福 社 
和 实现 区 域 可 持续 性 具有 重要 的 意义 。 

日 前 ,国内 外 学 者 对 生物 多 样 性 、 产 水 和 固 矶 每 生态 系统 服务 之 间 的 权衡 和 协同 关系 做 了 一 定 的 探索 。 
Hall 等 “利用 GIS 空间 分 析 方法 探讨 了 增加 林木 覆盖 面积 对 森林 碳 储量 和 生物 多 样 性 的 影响 ,并 且 固 碳 与 生 
物 多 样 性 存在 着 复杂 的 关系 ; Hatanaka 等 “ 通过 对 树龄 分 级 ,分 析 了 固 碳 与 生物 多 样 性 之 间 的 权衡 关系 ; 
Bradshaw 等 ”从 碳 排放 价格 方面 研究 了 澳大利亚 不 同 碳 排放 价格 对 生物 多 样 性 的 影响 ;Rodrfguez 等 ”指出 
造林 虽然 可 以 增加 固 矶 量 ,调节 气候 ,但 也 会 增 大 共 散 量 ,减少 径流 量 ; 在 实践 中 需要 对 碳 价 格 和 水 价格 做 出 
权衡 ,从 而 确定 造林 面积 ;Onaindia 等 ”在 不 同 碳 价格 基础 上 探究 了 生物 多 样 性 “ 固 碳 和 水 供应 之 间 的 权衡 关 
系 ,提出 在 经 济 的 趋势 下 增加 人 工 林 会 减少 水 的 提供 并 危害 生物 多 样 性 。 李 晶 等 " 利用 相关 系数 和 空间 制 
图 研究 了 天天 经 济 区 的 NPP . 固 碳 释 氧 \, 水 文 调 水 土 保 持 等 多 种 服务 的 权衡 协同 关系 ,并 指出 调节 型 生态 
系统 服务 与 供给 型 生态 系统 服务 之 间 存 在 权衡 天 系 然而 ,国内 外 对 生态 系统 服务 之 间 权 衡 与 协同 关系 的 定 
量化 及 其 驱动 机 制 研究 还 远 远 不 够 ,多 重生 态 系 统 服 务 之 间 非 线性 动态 关系 仍然 非 稼 模糊 。 此 外 , 随 春 西部 
经 济 大 开发 和 “一 市 一 路 ”等 经 济 战略 的 实施 ,渭河 流域 城市 化 进程 加 快 , 人 们 对 生态 系统 服务 的 需求 也 急剧 
增加 ,导致 了 严重 的 生态 环境 问题 ,具体 表现 为 水 资源 短缺 .环境 污染 ,物种 多 样 性 减少 .生境 破碎 等 "” 。 
因此 ,加 强 态 系统 服务 之 间 及 其 与 人 类 发 展 和 生态 环境 的 权衡 分 析 , 为 区 域 经 济 发 展 和 生态 建设 提供 决策 依 
据 是 非 第 有 必要 的 。 

鉴于 未 来 的 不 确定 性 与 现实 决策 的 需要 ,土地 利用 情景 预测 是 目前 国际 上 进行 生态 系统 服务 权衡 与 协同 
分 析 的 重要 方法 “而 有 关 滑 河流 域 的 研究 中 ,对 于 未 来 土地 利用 变化 情景 的 设置 多 以 土地 利用 现状 为 
基础 , 何 单 设置 不 同 土地 利用 类 型 面积 的 增 减 变化 , 较 少 结合 研究 区 域 的 实际 和 发 展 规 划 。 因 此 ,本 文 以 渭河 
流域 (天天 段 ) 为 研究 对 象 ,通过 设 定 经 济 发 展 优先 环境 保护 优先 和 两 者 兼顾 等 一 系列 土地 利用 情景 ,定量 
模拟 了 未 来 多 种 生态 系统 服务 的 质量 与 空间 格局 变化 ,并 进一步 通过 权衡 与 协同 分 析 , 结 合 当 地 的 发 展 目标 
或 利益 需求 ,得 出 最 优 的 土地 利用 类 型 分 布 ,为 该 区 域 的 可 持续 发 展 和 生态 文明 建设 提供 决策 文 持 。 


1 研究 区 概况 和 数据 来 源 


1 1 研究 区 概况 

渭河 是 黄河 的 最 大 的 文 流 ,主要 流 经 甘肃 天 水 .陕西 宝鸡 .咸阳 西安 等 地 ,在 陕西 潼关 县 汇 人 黄河 ,全 长 
818 km ,多 年 平均 径流 量 75.7 xl10 m- 。 涓 河流 域 地 形 特点 为 西高 东 低 ,可 分 为 东西 二 部 : 西 为 黄土 丘陵 沟 餐 
区 , 东 为 关中 平原 区 ; 南 有 东西 走向 的 秦岭 横 百 , 北 有 六 盘山 屏障 。 渭 河 径流 年 际 变 化 表现 为 南部 小 ,北部 大 ; 
季节 变化 明显 ,干流 以 秋季 流量 最 大 , 约 占 年 径流 的 38% 一 40%。 渭 河流 域 关 天 有 段 涵盖 6 市 1 区 ,从 西向 东 分 
别 为 :天 水 市 .宝鸡 市 ,杨凌 示 放 区、 西安 市 .咸阳 市 .铜川 市 .渭南 市 .商洛 市 部 分 区 县 (图 1)。 此 区 域 为 陕西 
省 的 行政 中 心 和 经 济 枢纽 ,人 口 密度 和 城市 化 程度 较 高 。 消 河 承 担 厦 经济 区 的 生产 、 生 活用 水 以 及 生态 补水 
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的 水 域 功 能 。 


图 例 
e 城市 e 采样 点 河流 [| 市 级 
m 子 流域 [ 渭河 流域 (关中 天 水 经 济 区 段 ) 


105° 106° 107° 108° 109° 110°E 


1 渭河 流域 关 天 段 区 域 图 


Fig.1 The location of Cuantian section in Weihe River basin 


1.2 ”数据 与 处 理 

本 研究 选用 的 数据 主要 包括 渭河 流域 关 天 段 的 基础 地 理 信息 数据 (行政 区 县 .河流 .道路 等 ) .遥感 影像 
图 DEM 数据 .土地 利用 数据 (2000 年 .2005 年 .2010 年 2013 年 ) 气象 数据 (降水 温度、 风速 .相对 湿度 等 )， 
统计 数据 ( 人口 .GDP 等 ) 。 野 外 采样 与 实地 调查 数据 ,沿路 线 西 安 -商洛 -铜川 -渭南 -宝鸡 -天 水 -杨凌 -咸阳 - 西 
安 环 绕 渭 河流 域 关 天 段 ,共和 采集 90 个 土壤 样 点 数据 (图 1) 。 通 过 重 铬 酸 钾 氧 化 -外 加 热 法 测算 采样 点 的 土壤 
有 机 碳 含 量 “ ,并 与 碳 循环 模型 反 演 土壤 呼吸 计算 得 出 的 土壤 有 机 碳 含量 进行 对 比 验 证 ,相关 系数 达到 
0.7391。 并 通过 实地 调查 及 谷歌 地 球 影像 验证 遥感 解 译 的 土地 利用 类 型 图 并 修改 与 实际 不 符合 的 部 分 。 


2 研究 方法 


2.1 固 碳 估算 

本 文 研究 的 固 碳 主 要 分 为 地 上 与 地 下 两 部 分 ,地 上 部 分 由 净 初 级 生产 力 ( NPP ) 推算 得 出 ,每 生死 干 物质 
会 固定 1.63 kg 的 -C0;; 地 下 部 分 主要 关注 土壤 有 机 碳 。 

(1) 地 上 部 分 采用 -CASA( Carnegie-Ames-Stanford Approach ) 模型 计算 植被 净 初 级 生产 力 , 进 而 计算 出 地 
上 固 碳 量 o 兆 初 级 生产 力 的 估算 由 植被 吸收 的 光合 有 效 辐射 ( APAR) 和 光 能 转化 率 (e) 两 个 变量 确定 。 

NPP(«,t) = APAR(x,t) x e(x,t) (1) 

AH sac RREI; RRN; APAR (x,t) 表示 空间 位 置 * 在 1 月 吸收 的 光合 有 效 辐射 (MJ m2") ;e 
(x,t) DRT 在 1 月 的 光 能 转化 率 (gC/MJ)。 

(2) 地 下 部 分 选用 了 由 周涛 等 "7 改进 的 碳 循环 过 程 异型 来 反 演 土壤 基础 呼吸 ,进而 与 全 国 第 二 次 土壤 普 
查 的 士 壤 有 机 碳 数据 进行 回归 分 析 , 得 到 土壤 呼吸 与 土壤 有 机 碳 含量 的 回归 方程 。 碳 循环 过 程 模型 的 计算 公 
式 如 下 : 


NPP 
式 中 ,4; 为 土壤 基础 呼吸 ;NPP 为 净 初 级 生产 力 ;b 为 温度 敏感 常数 因子 。 
2.2 生物 多 样 性 估算 
本 文中 应 用 InVEST 软件 的 生物 多 样 性 模块 测算 生境 质量 指数 评价 生境 质量 。 计 算 结果 范围 为 :0 一 1， 


(2) 
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值 越 大 代表 栖息 地 质量 越 蜗 ,生物 多 样 性 质量 越 品 。 计 算 公式 为 

Q,, = A — (Da De tk))) (3) 
式 中 ,0, 是 土地 利用 类 型 j 中 栅 格 x 的 生境 质量 ;D, 是 土地 利用 类 型 j 栅 格 x 的 生境 胁迫 水 平 ;k 是 缩放 参数 
(HRO) o 


R Y R d,, 
Dad a es (4) 


max 


IF d 为 顶 格 * 与 y 之 间 的 线性 距离 ;qd,,,, 为 胁迫 因子 最 大 有 效 威胁 距离 5 WS AF AE SB NA x 
的 可 达 性 水 平 ;5; 为 土地 利用 类 型 7 对 胁迫 因子 的 敏感 性 。 
2.3 ” 产 水 估算 

本 文 构建 了 渭河 流域 的 SWAT 模型 ,并 根据 研究 区 面积 .土地 利用 类 型 .土壤 属性 等 实际 情况 , 设 定 最 小 
集 水 面积 靖 值 为 13000 hm? ,将 渭河 流域 关 天 段 划分 为 157 个 子 流域 (图 1)。 产 水 服务 是 基于 SWAT 模型 水 
文 模块 测算 的 ,其 计算 公式 如 下 : 

WYLD =SURQ+LATQ+GWQ-TLOSS-PA (5) 
式 中 ,WYLD 为 总 产 水 量 ,是 指 时 间 步 长 内 进入 主 河道 的 总 水 量 (mm) ;SURQ 为 时 间 步 长 内 地 表 径 流 对 主 河 
道 总 径流 的 贡献 量 (mm) ;LATQ 为 时 间 步 长 内 , 侧 向 流 对 河川 径流 的 贡献 量 (mm) ;GWQ 为 时 间 步 长 内 ,地 下 
径流 对 主 河道 总 径流 的 贡献 量 (mm) ;TLOSS 为 河床 传输 的 水 损失 量 (mm) ;Ph4 为 池塘 截留 量 (mm)。 

为 了 验证 构建 模型 的 准确 性 ,研究 选择 1990 年 至 2013 年 咸阳 水 文 站 的 实测 径流 数据 对 SWAT 模型 进行 
参数 的 率 定 和 验证 。 其 中 1990—1999 年 为 率 定期 /以 2000 年 .2005 年 .2010 年 .2013 年 的 径流 数据 进行 验 
证 。1990 一 1999 年 逐 月 实测 径流 量 与 模拟 值 分 布 接近 ,相关 性 比较 好 ,相关 系数 尺 为 0.78,NES 值 为 0.71 , 均 
KF 0.5; 验 证 阶段 2000 年 至 2013 年 ,这 4 年 的 模拟 值 和 实测 值 吻合 度 也 比较 好 ,相关 系数 尺 为 0.73( 图 2)。 
这 说 明 构 建 的 SWAT 模型 在 渭河 流域 具有 良好 的 适用 性 ,径流 模拟 精度 很 高 ,可 以 用 来 做 进一步 的 研究 工作 。 
2.4 土地 利用 情景 预测 

根据 流域 历史 土地 利用 转移 规律 ,结合 未 来 发 展 规 划 , 构 建 了 流域 未 来 3 种 土地 利用 发 展 情景 ,利用 
IDRISI 软件 的 CA-Markov 模块 ”对 流域 2050 年 的 土地 利用 进行 模拟 。 不 同情 景 规划 下 ,对 高 程 、 降 水 .人 
O 道路、 水系 等 驱动 及 限制 因子 设置 不 同 参数 与 权重 值 ,得 到 2050 年 土地 利用 类 型 图 (图 3) 。 

情景 1 为 计划 情景 ,土地 利用 变化 按 原 有 趋势 持续 演变 ,并 参考 渭河 流域 土地 利用 总 体 规划 ,调整 未 来 土 
地 利用 类 型 面积 和 分 布 s 

情景 2 为 开发 情景 ,在 现 有 土地 利用 的 基础 上 推进 城市 化 建设 ,城镇 边缘 的 耕地 , 稀 玻 洪 木 林地 大 量 转 化 
城镇 建设 用 地 。 

情景 3 为 保护 情景 ,限制 城镇 化 发 展 , 生 态 环境 保护 为 第 一 目的 。 严 格 限制 城镇 建设 .工业 用 地 侵占 耕地 
和 水 域 ,对 坡度 大 于 20?" 的 耕地 进行 退耕 还 林 还 草 措施 。 

2.5 .生态 系统 服务 优化 

在 本 文中 ,研究 引入 生产 可 能 性 边界 方法 对 未 来 生态 系统 服务 的 进行 优化 ,并 得 到 最 优生 态 系统 服务 对 
应 的 十 地 利用 类 型 。 生产 可 能 性 边界 ; 也 称 为 由 累 托 曲线 ( Pareto frontiers ) , 是 指 在 资源 量 固定 不 变 的 情景 
下 ,所 能 够 产 出 的 一 种 或 多 种 产品 的 量 的 总 和 - 。 数 学 方式 表达 出 来 即 为 一 条 坐标 轴 中 的 曲线 ,可 以 用 来 反 
映 潜 力 与 过 度 问 题 。 如 果 产 品 组 合 点 在 曲线 的 下 方 ,说 明 资 源 利 用 不 充分 ,存在 优化 潜力 ; 奋 组 合 点 落 在 曲线 
上 , 则 为 最 优 配置 ; 知 落 在 曲线 以 外 ,那么 是 无 法 达到 的 。 下 面 以 生物 多 样 性 与 产 水 两 种 生态 系统 服务 为 例 ， 
说 明 制 作 两 者 之 间 的 由 累 托 效率 曲线 的 方法 。 

首先 ,对 生物 多 样 性 栅 格 图 层 与 产 水 图 层 相 除 得 到 比值 图 层 ,比值 图 层 中 每 个 单元 格 代 表 对 应 地 理 位 置 
的 生物 多 样 性 与 产 水 的 比值 ;然后 ,按照 生物 多 样 性 / 产 水 的 大 小 对 比值 图 层 每 个 单元 格 做 升序 排列 ;最 后 , 按 
照排 列 的 升序 顺序 ,依次 分 别 对 单元 格 对 应 地 理 位 置 的 生物 多 样 性 质量 和 产 水 量 进行 累计 求 和 ,并 按照 结 
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Fig.2 Simulated and realistic monthly runoff data in calibration and validation period 
绘制 曲线 , 即 得 到 生物 多 样 性 - 产 水 之 间 的 帕 累 托 效 率 曲 线 。 
3 结果 与 分 析 


3.1 生态 系统 服务 时 空 变化 

运用 相关 模型 ,测算 并 制图 渭河 流域 2000—2013 年 的 产 水 , 固 左 和 生物 多 样 性 (图 4)。 渭 河流 域 产 水 量 
呈现 逐年 增多 的 趋势 。2000 年 所 有 子 流域 的 平均 产 水 量 为 110.847 mm, 2005 年 平均 产 水 量 比 2000 年 小 幅度 
增加 了 31.5 mm, 到 子 2010 年 和 2013 年 平均 产 水 量 大 幅 增 长 ,分 别 是 247.91 mm 和 275.42 mm。 主 要 是 因为 
2010 年 和 2013 年 的 降水 量 增加 ,尤其 是 在 7、8 月 份 的 高 强度 降雨 ,造成 河川 径流 量 也 大 幅 增 长 。 在 空间 上 
分 析 ; 除 了 2000 年 的 产 水 分 布 是 下 游 最 高 ,中游 略 小 于 上 游 外 ,其 他 几 个 年 份 都 是 产 水 量 下 游 远 远 高 于 中 游 
和 和 上游; 二 游 的 产 水 量 最 小 。 除 此 以 外 ,4 个 年 份 的 最 高 产 水 量 均 出 现在 下 游 的 2 号 子 流域 。 通 过 测算 地 上 、 
地 下 固 碳 量 , 发 现 2000 年 至 2013 年 固 碳 量 空间 分 布 特征 总 体 变 化 不 大 ,平均 固 碳 量 最 高 的 区 域 在 62 号 子 流 
域 ,最 低 值 在 3 号 子 流域 ,其 中 2000 年 是 最 低 值 为 27.6 Vhm2。62 号 子 流域 及 其 相 邻 区 域 拥 有 较 高 的 平均 固 
碳 量 是 因为 它们 的 土地 利用 类 型 基本 都 是 林地 ,植被 覆盖 率 高 的 地 方 固 碳 量 会 比较 高 。 关 中 平原 地 区 子 流 域 
国 碳 量 相对 较 少 ,秦岭 地 区 子 流 域 固 破 量 较 大 ;局 部 变化 比较 明显 ,个 别 子 流域 从 2000 年 至 2013 年 的 固 碳 量 
有 明显 变化 ,降水 量变 化 及 子 流 域内 土地 利用 类 型 变更 可 能 是 导致 固 碳 量变 化 的 主要 原因 。2000 年 至 2013 
年 渭河 流域 关 天 段 生物 多 样 性 质量 呈 下 降 趋势 。 渭 河流 域 关 天 段 的 生物 多 样 性 得 分 范围 为 0.29 一 0.73 ,城市 
建设 用 地 的 得 分 普遍 较 低 ,秦岭 山地 的 得 分 较 高 ,这 表明 秦岭 地 区 的 生态 退化 程度 低 ,生物 多 样 性 质量 高 。 如 
位 于 秦岭 山地 的 157 号 子 流域 ,生态 威胁 因子 数量 较 少 ,生境 与 威胁 因子 距离 较 远 ,生物 多 样 性 质量 高 。 而 关 
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Fig.3 The land use patterns for 2000, 2010 and 2050 in Guantian section of Weihe River basin 


中 平原 生境 质量 得 分 在 0.38 左右 ,该 区 域 受 人 为 影响 较 大 , SELL Se CE So Bae SE AB E PE AE AS WES PED, 
生境 质量 下 降 。 
3.2 ”生态 系统 服务 之 间 相 关 性 

本 文 以 子 流域 为 研究 对 象 ,运用 R 语言 编程 实现 对 子 流域 样 点 的 相关 性 分 析 '"… ALS 的 各 图 左下 角 为 
散 点 图 ,右上 角 为 用 于 展示 相关 系数 的 饼 图 ; 散 点 图 与 饼 状 图 共同 表达 生态 系统 服务 之 间 相 关 性 的 情况 。 饼 
状 图 填充 大 小 代表 相关 系数 值 , 填 充满 表 示 相 关系 数 为 1。 颜色 深浅 代表 值 的 大 小 , 值 越 大 颜色 越 深 , 蓝 色 代 
表 正 相关 ,粉色 代表 负 相 关 。 硅 为 正 相关 关系 则 顺 时 针 填 充 饼 图 , 奉 为 负 相 关 关 系 则 逆 时 针 填 充 饼 图 。 由 图 
5 可 以 看 出 ,2000 年 至 2013 年 , 固 砚 与 生物 多 样 性 均 表 现 为 正 相 关 ,2010 年 相关 系数 最 大 , 约 为 0.78 ,2013 年 
最 小 , 约 为 0.63。 固 碳 与 产 水 均 表 现 为 负 相关 ,2005 年 相关 系数 最 大 , 约 为 0.45 ,2000 年 最 小 , 约 为 0.13; 生 物 
多 样 性 与 产 水 之 间 为 外 相关 ,2010 年 相关 系数 最 大 , 约 为 0.26 ,2000 年 最 小 , 约 为 0.08。 

运用 Python 语言 绘制 固 碳 产 水 和 生物 多 样 性 三 者 之 间 的 PPF 曲线 ,进而 对 渭河 流域 关 天 上 段 这 一 固定 资 
源 量 的 生态 系统 进行 生态 系统 服务 之 间 权 衡 与 协同 关系 的 定量 研究 。 如 图 6 所 示 , 固 碳 与 生物 多 样 性 之 
间 呈 明显 协同 关系 ,从 点 a >b 和 c 一 d, 固 碳 增加 量 都 是 5S0x10"t, 而 生物 多 样 性 分 别 增 加 了 0.63x10" 和 0.72 
x10" ;这 说 明 固 碳 量 累计 越 多 ,对 生物 多 样 性 质量 越 有 利 。 图 6 所 示 , 固 碳 与 产 水 以 及 生物 多 样 性 与 产 水 之 间 
的 PPF 曲 线 表 现 为 回 外 “点 ”的 曲线 ,存在 此 消 彼 长 的 权衡 关系 。 以 图 6b P a.b cd 4 个 点 进行 固 碳 与 产 水 
之 间 的 权衡 关系 分 析 。 从 a 点 到 b 点 ,变化 幅度 比较 小 ,平均 和 斜率 为 -204.7/50, 即 增加 50x10° +t 的 固 碳 量 
会 减少 204.7x105t 的 产 水 量 ;从 c 点 到 d 点 的 斜率 为 -541.1/50, 即 增加 50x10st 的 固 碳 量 会 减少 341.1x105t 
的 产 水 量 。 研 究 称 固 碳 为 产 水 的 机 会 成 本 ,不同 阶段 减少 产 水 量 所 付出 的 计划 成 本 是 有 差异 的 ,从 a b 和 
>d 付出 相同 的 成 本 带 来 的 效益 是 不 同 的 ,显然 从 cd 获得 的 收益 比较 大 ,这 表明 固 碳 量 积 累 越 多 越 有 利于 
保水 。 产 水 与 生物 多 样 性 之 间 的 由 累 托 效率 曲线 的 斜率 同样 也 是 从 右 向 左 逐 渐变 大 。 
3.3 生态 系统 服务 优化 配置 

本 文 模拟 了 2050 年 的 土地 利用 类 型 情景 ,并 测算 了 3 种 情景 下 的 固 碳 、 生 物 多 样 性 和 产 水 3 种 生态 系统 
服务 量 。 为 了 使 未 来 的 3 种 生态 系统 服务 能 够 达到 最 大 化 ,选取 2050 年 3 种 情景 下 的 固 碳 、 生 物 多样 性 与 产 
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Fig.4 The spatial distribution of water yield, carbon sequestration and biodiversity from 2000 to 2013 


水 3 种 生态 系统 服务 的 最 大 值 绘制 最 优 帕 累 托 效率 曲线 ,代表 未 来 土地 资源 所 能 提供 的 最 优 的 生态 系统 服务 
总 量 。 并 将 2013 年 的 土地 利用 资源 配置 下 的 生态 系统 服务 组 合 情 况 展 点 到 坐标 系 中 ,发 现 2013 年 生态 系统 
服务 的 组 合 点 均 在 最 优 帕 累 托 效率 曲线 下 方 ,说 明 资 源 配置 存在 优化 的 潜力 。 下 面 我 们 将 以 固 碳 与 生物 多 样 
性 为 例 , 探 究 这 两 种 生态 系统 服务 最 优 帕 累 托 效率 曲线 下 的 土地 利用 优化 。 

在 图 7 中 ,点 b(150;1:88) 为 2013 年 的 固 克 和 生物 多 样 性 ,所 对 应 的 土地 利用 类 型 在 空间 分 布 上 包括 渭 
河流 域 关 天 段 的 上 游 和 关中 平原 的 大 部 分 区 域 ,耕地 面积 为 2.27 «10° hm ,林地 面积 为 5.19x10” hm ,草地 面 
积 为 9.51x10 hm ,水域 面积 为 1.02x10” hm? ,城镇 用 地 面积 为 9.13x10” hm ”。 最 优 帕 累 托 效率 曲线 上 的 点 c 
(150,2.16) 和 点 d(132,1.88) 是 2050 年 最 优 配置 下 的 固 兢 和 生物 多 样 性 。 对 于 存在 优化 潜力 的 点 b(150， 
1.88) 可 以 采取 两 种 极端 的 方式 到 达 最 优 帕 累 托 效率 曲线 ,一 种 是 通过 上 移 至 点 ec(150,2.16) , 固 碳 量 保 持 不 
变 , 增 加 生物 多 样 性 ;需要 将 b 点 对 应 的 土地 利用 优化 成 ce 点 所 对 应 的 土地 利用 类 型 及 分 布 。e 点 对 应 的 土 
地 利用 分 类 中 ,耕地 面积 为 2.45x10° hm? ,林地 为 5.76x105 hm? ,草地 为 9.34x105 hm2 ,水 域 为 1.05x10* hm’, 
城镇 用 地 为 1.13x10 hm, XTE b 点 地 类 面积 ,e 点 中 只 有 草地 的 面积 减少 ,其 余地 类 的 面积 均 有 所 增加 。b、 
c 两 者 之 间 的 空间 分 布 差异 较 大 ,b 点 中 南部 的 地 类 去 除 ,而 东南 与 西南 方向 增加 了 区 域 。 另 一 种 方式 是 右 移 
至 点 d(132,1.88) ,生物 多 样 性 保持 不 变 , 减 少 固 碳 量 ;需要 将 所 对 应 的 土地 利用 转变 成 d 点 对 应 的 土地 利用 
类 型 及 分 布 。d 点 对 应 的 耕地 面积 为 2.25x10* hm ,林地 为 4.32x10” hm? ,草地 为 8.14x10” hm ,水域 为 1.01x 
10 hm? ,城镇 用 地 为 1.03x10” hm? 。 实 际 中 ,只 要 能 够 通过 改变 土地 利用 方式 使 更 加 接近 这 条 曲线 , 即 为 实现 
了 生态 系统 服务 的 优化 。 

当 2013 年 的 累计 固 碳 量 达到 最 大 值 时 ,就 是 坐标 系 中 的 点 a, 此 时 所 对 应 的 土地 利用 分 布 就 是 整个 研究 
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Fig.5 The correlations of ecosystem services in Cuantian section of Weihe River basin from 2000 to 2013 


区 渭河 流域 ( 关 天 段 ) 2 研究 可 以 根据 曲线 对 整个 研究 区 的 土地 利用 类 型 进行 优化 ,使 生态 系统 服务 到 达 最 
优 值 , 即 达到 曲线 上 的 点 A。a\A 点 所 对 应 下 的 土地 利用 类 型 及 空间 分 布 情况 如 图 8 所 示 。 综 上 ,我 们 可 以 
利用 最 优 幅 素 托 效 率 曲 线 ,权衡 不 同 的 目标 和 生态 服务 需求 ,进行 定量 化 的 土地 利用 优化 ,从 而 获取 最 大 化 的 
生态 系统 服务 价值 。 


本 文 引 入 生产 可 能 性 边界 对 固 左 、 生 物 多 样 性 和 产 水 之 间 协 同 与 权衡 关系 进行 了 定量 化 地 动态 分 析 , 并 
利用 R 和 Python 语言 对 三 者 之 间 的 协同 与 权衡 关系 进行 了 更 加 直观 清晰 地 表达 ,而 不 是 仅仅 停留 在 相关 性 
等 宏观 分 析 。 从 科学 研究 层面 讲 , 多 重生 态 系统 服务 之 间 及 其 与 地 理 环 境 之 间 复 杂 关 系 的 变化 特征 和 作用 机 
理 还 需 进一步 的 归纳 和 总 结 。 在 现 有 工作 的 基础 上 ,进一步 研究 生态 系统 服务 协同 与 权衡 关系 动态 变化 的 驱 
动机 制 和 探讨 地 理 环境 空间 异 质 性 下 的 尺度 效应 是 本 文 下 一 步 的 研究 方向 ;从 实际 决策 层面 讲 , 人 类 活动 
(如 城市 化 ,退耕 还 林 等 ) 会 导致 地 表 类 型 发 生 剧 烈 改 变 , 而 土地 利用 /罗兰 变化 是 引起 生态 系统 服务 及 其 彼 
此 之 间 的 权衡 关系 发 生变 化 的 主要 驱动 因素 ”” 。 例 如 ,本 文 的 研究 结果 表明 , 产 水 与 固 碳 之 间 存 在 权衡 关 
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6 ”渭河 流域 关 天 上 段 固 碳 、 生 物 多 样 性 和 产 水 之 间 的 生产 可 能 性 边界 (PPF) 曲线 


Fig.6 The trade-off and synergies between carbon sequestration, biodiversity and water yield in Guantian section of Weihe River basin 


系 ,相关 性 系数 在 2000 年 为 0.13, 但 在 2013 年 相关 系数 出 现 显 著 增 长 , 约 为 0.45, 这 主要 是 归 因 于 降水 量 的 
增加 和 林地 的 减少 。2005 年 与 2000 年 相 比 ,土地 利用 类 型 也 发 生 了 明显 变化 ,林地 面积 减少 62214 hm , 城 
镇 用 地 增加 了 18504‘hm ,植被 减少 直接 影响 碳 固 定 , 同 时 城市 化 也 加 快 了 矶 排放 和 产 流 的 速率 。2005 年 与 
2000 年 相 比 ,降水 量 增加 了 了 约 70 mm, 这 直接 导致 产 水 量 增加 。 在 经 济 迅 速 发 展 的 同时 ,人 类 活动 对 自然 环 
境 造 成 了 巨大 的 负面 影响 ,为 了 实现 永 续 发 展 ,必须 在 经 济 发 展 和 生态 保护 之 间 找 到 平衡 点 。 这 就 在 客观 上 
要 求 强化 对 多 种 服务 之 间 及 其 与 生态 环境 之 间 的 权衡 分 析 , 从 而 为 区 域 发 展 和 生态 建设 提供 决策 依据 ;从 具 
体 应 用 层面 讲 , 生 态 系 统 服 务 的 情景 模拟 与 权衡 分 析 是 针对 性 解决 生态 系统 和 土地 利用 管理 中 存在 问题 的 有 
效 手段 ;可 以 为 区 域 国 土 规划 生物 多 样 性 保护 和 生态 补偿 等 提供 科学 依据 “i 。 例 如 ,退耕 还 林 或 重 造 森 
林 会 丰富 生物 多 样 性 ,增强 气候 调节 和 和 森林 水 质 净化 等 功能 ,但 同时 可 能 会 增加 蒸 散 发 ,减少 产 水 量 ””” ,并 
影响 农业 生产 活动 中 食物 和 原材料 的 供给 等 ”。 本 人 研究 结果 也 表明 ,在 流域 的 实际 管理 过 程 中 ,利益 相关 考 
可 以 通过 优化 土地 利用 使 固 兢 量 和 生物 多 样 性 同时 增加 ,或 者 通过 权衡 需求 使 固 碳 量 或 者 产 水 量 其 中 一 个 增 
加 , 男 一 个 减少 。 

本 人 研究 中 设 定 经 济 发 展 优先 、 环 境 保 护 优先 和 两 者 莱 顾 一 系列 土地 利用 情景 ,可 以 定量 模拟 未 来 多 种 生 
态 系统 服务 的 质量 与 空间 格局 变化 ,进一步 通过 权衡 与 协同 分 析 ,结合 当地 的 发 展 目标 或 利益 需求 ,选择 最 优 
方案 并 提出 相应 的 优化 管理 措施 。 但 未 来 土地 利用 的 发 展 受 多 方面 的 影响 和 制约 ,并 存在 一 些 不 可 控 因 素 
(例如 政策 导向 因素 ) ,可 能 需要 模拟 几 十 甚至 上 百 种 情景 来 还 原 现实 。 而 情景 只 是 对 可 能 的 未 来 的 描述 , 鉴 
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7 未 来 情景 下 生态 系统 服务 之 间 的 权衡 与 协同 关系 及 土地 利用 优化 


Fig.7 Trade-off and synergy between ecosystem services in future scenarios 
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8 最 优 帕 累 托 效率 曲线 对 应 的 土地 利用 分 布 图 


Fig.8 The land use patterns under the optimal Pareto frontier 
于 未 来 的 不 确定 性 ,无疑 会 影响 未 来 生态 系统 服务 的 模拟 精度 ,增加 权衡 与 协同 分 析 的 不 确定 性 … 。 


此 ,在 接 下 来 的 研究 中 进一步 完善 土地 利用 情景 设置 ,加强 对 不 确定 性 的 定量 化 分 析 及 空间 制图 化 表达 ,对 于 
生态 系统 管理 和 土地 利用 决策 具有 重要 的 意义 。 


ait 


本 文 以 渭河 流域 天 天 上段 为 研究 对 象 ,在 子 流域 尺度 估算 了 2000 年 ,2005 年 ,2010 4F , 2013 年 的 生物 多 样 
性 、 固 矶 和 产 水 3 种 生态 系统 服务 ,并 对 生态 系统 服务 彼此 之 间 权 衡 和 协同 关系 进行 定量 化 的 动态 分 析 ; 并 通 
过 土地 利用 情景 预测 和 最 优 帕 累 托 效 率 曲线 优化 了 2050 年 最 优生 态 系统 服务 下 的 土地 利用 格局 。 结 果 表 
明 :(1)2000 一 2013 年 以 来 ,渭河 流域 关 天 上 段 总 产 水 量 呈 增加 趋势 ,主要 分 布 于 渭河 流域 下 游 地 区 ;秦岭 与 六 
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碳 能 力 差 , 受 人 类 活动 扰动 因素 较 多 ,生物 多 样 性 质量 低 。(2) 渭 河流 域 关 天 段 的 生态 系统 服务 呈现 权衡 加 
强 协同 略 减 的 趋势 。 具 体 而 言 : 固 兢 和 生物 多 样 性 之 间 为 协同 关系 , 且 14 年 来 相关 性 呈现 上 下 波动 ,总 体 上 
略微 减弱 的 趋势 ; 产 水 与 固 碳 为 权衡 关系 ,相关 性 在 2005 年 出 现 显 著 增 强 , 之 后 逐渐 减弱 ; 产 水 与 生物 多 样 性 
为 权衡 关系 ,相关 性 在 2000 一 2010 年 呈现 逐渐 增强 ,之 后 在 2015 年 略 有 减弱 。(3) 绘 制 了 2050 年 的 生态 系 
统 服 务 下 的 最 优 帕 累 托 效率 曲线 并 得 出 对 应 的 土地 利用 优化 图 ,通过 转变 土地 利用 类 型 使 生态 服务 功能 达到 
规划 目标 。 这 些 有 关 生 态 系 统 服务 协同 与 权衡 关系 及 其 未 来 最 大 生态 系统 服务 指导 下 的 土地 利用 优化 的 研 
究 结 果 , 为 延 河流 域 土 地 资源 管理 和 协调 可 持续 发 展 提供 了 必要 的 科学 依据 。 
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